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We experimentally study particle scale dynamics during segregation of a bidisperse mixture under 
oscillatory shear. Large and small particles show an underlying asymmetry that is dependent on 
the local particle concentration, with small particles segregating faster in regions of many large 
particles and large particles segregating slower in regions of many small particles. We quantify the 
asymmetry on bulk and particle scales, and capture it theoretically. This gives new physical insight 
into segregation and reveals a similarity with sedimentation, traffic flow and particle diffusion. 

PACS numbers: 45.70.Mg, 05.45.-a, 47.57.Gc 


The natural tendency of granular media to self- 
organize when agitated or sheared produces a rich diver¬ 
sity of complex and beautiful patterns [J-fll- Although 
it is counter-intuitive that the components of a heteroge¬ 
neous mixture will readily separate, this property has se¬ 
rious technical implications as the cause ofproduct non¬ 
uniformity in many industrial processes [^-|6| and also 
plays a pivotal role in the enhanced run-out of large scale 
geophysical granular flows, such as debris-flows, pyroclas¬ 
tic flows and snow avalanches EM- A firm knowledge of 
the segregation process is thus of universal importance. 

Although there has been considerable recent progress 
in developing continuum based segregation models for 
sheared granular flows EHH, the individual particle 
dynamics are still poorly understood. Discrete Particle 
Method (DPM) simulations 17-2Q1 produce a wealth of 
micro-scale information, but are models in themselves. It 
is vital to directly measure particle segregation dynam¬ 
ics in real experiments, but such an analysis is difficult 
with conventional techniques such as binning and side- 
wall observation [2ll - [24j . Non-intrusive imaging tech¬ 
niques, such as X-ray tomography 25] and refractive 
index-matched scanning (RIMS) |26l . |27| allow examina¬ 
tion of the interior of a granular medium, with RIMS 
recently developing into a powerful tool for examining 
monodisperse and bidisperse flows 29j-31|. In particu¬ 
lar, the work of Harrington et al. [31] on the emergence 
of granular segregation demonstrates how particle scale 
analysis can give new physical insights. 

In this Letter, we analyze particle scale dynamics dur¬ 
ing segregation of a bidisperse mixture under oscillatory 
shear. We find that the behavior of small and large par¬ 
ticles exhibits an asymmetry related to the local particle 
concentration, with small grains moving faster through 
regions of many large particles and large particles ris¬ 
ing slower through regions of many small particles. This 
asymmetry is quantified on both particle and bulk length 
scales, and it is shown how to incorporate the behavior 
within the theoretical framework. 

Methods.— A shearbox 51 mm deep and 37 mm wide 
is filled to a height h = 87 ± 3 mm with a bidisperse 


mixture of borosilicate glass spheres ( p v = 2.23 g/cm 3 ) 
with diameters di = 8 mm and d s = 4 mm. The larger 
particles are placed at the bottom, the surface flattened 
and the smaller particles placed on top. The sidewalls 
oscillate whilst remaini ng p arallel, applying a periodic 
shear q(t) = 70 sin(uft) [32| as shown in Figure [I| The 
corresponding shear rate 7 (t) = 7owcos(wi), frequency 
w = 2n/T rads -1 , period T = 13 s and strain ampli¬ 
tude 70 = tan9 max . The sidewalls displace to a maxi¬ 
mum angle 9 max = ±30°, giving a maximum shear rate 
of 70 = 7o^ and a maximum grain displacement ampli¬ 
tude A = hjo- The angle is decreased to 9 ma x = ±10° 
for the particle trajectory data in order to slow down the 
segregation and increase the temporal resolution. Non- 
dimensional time t = t/T corresponds to the number of 
elapsed cycles. We follow a sample using RIMS, with 
the index-matched liquid a mixture of benzylalcohol and 
ethanol (viscosity fi = 3 mPas) containing a fluorescent 
dye (rhodamine). The low viscosity of the interstitial liq¬ 
uid means that fluid drag forces are small compared to 
both gravitational forces and the applied shear (Stokes 
number > 1 [33]). The mixture is lit with a 532 nm laser 
sheet perpendicular to the oscillating walls, giving a stack 
of vertical cross sections. A scan is performed after each 
full oscillation with the shearbox in the upright position. 
The images were processed using convolution [34j to give 



FIG. 1: (color online). Left: The experimental setup. A raw 
data image is shown and a cross-section of a reconstructed 
sample with 3 mm and 6 mm beads. Right: cross-sections at 
different times during an experiment. 
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three dimensional particle positions, which are coarse¬ 
grained in order to determine a continuous volume frac¬ 
tion [35J]. Some sidewall effects exist, with small particles 
preferentially located near the stationary vertical walls, 
but this does not affect the overall segregation. The hor¬ 
izontal particle motion is diffusive, hence the concentra¬ 
tion is spatially averaged to give a uniform concentration 
in the x-y plane. We observe no convection rolls [36J, al¬ 
though for \4>max\ > 45° geometrical squeezing was seen 
to cause convection. 

Results. —The typical behavior is shown in Fig. [I] The 
initial state with large particles on the bottom evolves to 
a final state with large particles on top, because small 
particles sink and large particles rise. Interestingly, some 
large particles remain below when all the others have 
reached the top. These particles are not stuck but rise 
at a slower rate than the ones that have reached the top 
before them. Although this has been inferred before, it 
has not yet been explained 0]. 

We define a segregation time t s as the time needed for 
the vertical centers of mass of the two species to reach a 
steady state, as shown in Fig.j^a). We record t s for mix¬ 
tures with varying global volume fraction of small parti¬ 
cles $(%) = V s /{yi+V s ), while keeping the total mixture 
volume constant. FigureE^b) shows that t s scales linearly 
with <f>, i.e. with more small particles in the mixture the 
segregation is slower (20]. Similar trends were observed 
over the entire range of angles that can be accessed in 
our setup. This behavior is usually given a two-part ex¬ 
planation: At low 4>, small particles move slower when 
there are more small particles [37]. At high $, it takes a 
longer time for large particles to travel to the top when 
the layer of small particles above them is thicker 0,13. 
In both explanations the behavior of the other species 
is ignored. So how do these explanations combine at an 
intermediate <f>? A clue is given by 0] which reported 
that for a 'F = 50% mixture the transition from the state 
with small particles on top to a mixed state was faster 




FIG. 2: (a) Time evolution of the vertical center of mass po¬ 
sition (4 X]"=i z i) f° r large and small particles in $ = 25% 
(black), 50% (dark gray) and 75% (light gray) mixtures. 
dmax = ±30°. (b) Segregation time t s as a function of <f>; solid 
line is a fit for the symmetric model with S r = 0.016, while 
the dashed line is a fit for asymmetric model with S r = 0.030 
and k = 0.89. 


than the subsequent transition from the mixed state to 
the final segregated state. This points to two separate 
processes that are likely to be related to the distinct be¬ 
havior of small and large particles. 

Particle dynamics .—We are thus motivated to study 
a single small particle segregating in a mixture of large 
particles and a single large particle segregating in a mix¬ 
ture of small particles, which we refer to as d* = 0 + % 
and <I> = 100~% mixtures respectively. The trajectories 
of the two particles, shown in Fig. [3](a), are quite differ¬ 
ent: (i) the large particle segregates roughly three times 
slower than the small particle; and (ii) the large particle 
rises smoothly at an almost constant speed, whereas the 
small particle shows a stepwise motion with steps of the 
order of di. This suggests that the small particle falls 
through gaps in the large particle matrix under gravity, 
typically traversing just a single layer. 

In order to more precisely understand the nature of 
these trajectories, we study the displacement after f cy¬ 
cles: Ar(r) = r{t + t) — r(f). The root mean square 
displacement (RMSD) a(f) = ^ (Ar 2 (f)) is plotted 
in Fig. [3|b). The dynamics are diffusive (logarithmic 
slope 1/4) for both particles at short timescales and su¬ 
per diffusive (logarithmic slope 1/2) at longer timescales. 
The crossover length scale between the diffusive and seg¬ 
regation (super-diffusive) regimes for the small particle is 
approximately di , which corresponds to the typical seg¬ 
regation step size of the small particle. The crossover 



FIG. 3: (color online) Individual particle dynamics for small 
particles (blue, gray) and large particles (black) with dmax = 
±10°. (a) Vertical trajectories of a small particle segregating 
in a <f> = 0+% mix; and a large particle segregating in a 
<f> = 10CF% mix. Inset: Particle movement in the horizontal 
plane, (b) RMSD <r(f) for different mixtures (see legend), 
with the solid line a fit of a s = \JDqt + w'jr' 2 at <E> = 50% 
(shifted for clarity). The dotted lines show the slopes 1/4 
and 1/2. Inset: w„($) for large (v = l) and small particles 
(u = s). (c)-(e) ? for single cycles in $ = 0 + %, 50% and 
100“% mixtures respectively. 
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length scale for the large particle is lower, roughly 0.2 di 
(0.4d s ), and is likely to be related to the scale of the re¬ 
arrangements of the surrounding small particles. To con¬ 
firm this, we measure the RMSD per cycle ? = (A r 2 }, 
as shown in Figs.0Kc) and0Ke). The typical value of <r for 
the single large particle lies just below e for the surround¬ 
ing small particles [Fig. 0Ke)]. Although the displace¬ 
ments g for the single small particle experiences large 
variations, as a result of falling through layers, the mean 
value is of the same order as that of the surrounding large 
particles [Fig.0Kc)]. 

The plot of <t(t) for a <f> = 50% mixture in Fig. 0Kb) 
shows that the curves lie between those for <f> = 0 + % 
and $ = 100“%, but with a comparable amount of 
segregation. Fitting each of the curves with = 

\/Dot + w 2 t 2 with diffusion coefficient Dq allows us to 
examine the segregation velocities w v for large (u = l) 
and small (v = s) particles at different <f>. The inset of 
Fig. 0Kb) shows that w s (<I>) decreases with increasing $, 
whereas wi( ( h) increases to a maximum at $ = 50% and 
then decreases, although not to zero, at $ = 100~%. To 
understand the peak in wi(&), we plot ? for <f> = 50% 
in Fig. [3jd) . The values of g for both small and large 
particles increase with respect to the <f> = 100“% mix, 
however, the large particle movement is still less com¬ 
pared to the small particles. 

At this point we can hypothesize an explanation for 
the trend in Fig. [2jb) : the individual dynamics of small 
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FIG. 4: (color online), (a)-(b) Conditional probabilities 

P(Azi\(j)) and P{Az s \ij>). White curves are (A zi) and (A z s ). 
(c)-(d) \(Azi)\/A and \(Az a )\/A as a function of <j>, with error- 
bars indicating the standard error of the mean. Dashed and 
solid lines are plots of Eq. © for quadratic and cubic flux 
functions F{<j >) with S r = 0.008 and S r = 0.015 respectively. 
The values of S r were scaled to account for the lower shear 
rate 70 CJ at 6 max = ± 10 °. 


and large grains have a different significance on the over¬ 
all segregation dynamics at different <f>. At high $, the 
significant dynamics are of the ‘slow’ large particle, which 
are governed by the scale of rearrangements of the sur¬ 
rounding small particles. At low $, it is the ‘fast’ small 
particle which is significant as it can make big segregation 
steps between large particle layers. At an intermediate $ 
both processes combine; small particles slow down, be¬ 
cause layering disappears, while large particles speed up, 
because the scale of rearrangements increases. 

Displacement statistics .—To study this behavior at the 
particle scale for each species {v = l,s ), we measure 
the conditional probabilities P(Az v \<f>) of the vertical dis¬ 
placement A z v given that the local small particle volume 
fraction is </>. Note that shear-gradients [38| do not play 
a role, because of the linear shear profile that is applied. 
Here, <fi = 0 corresponds to regions of only large parti¬ 
cles and <J) = 1 to only small particles. The results in 
Figs. UK a) and 01(b) demonstrate that large particles are 
less likely to segregate at high </> compared to small par¬ 
ticles segregating at low <f>. In the following, we will refer 
to this as “asymmetry”. Similar to the data for in 

the inset of Fig. 0Kb), we observe that the large particles 
have their greatest displacement at an intermediate <f>. 

The effect of asymmetry at a mesoscale can be seen in 
the temporal development of t) for a $ = 50% mix¬ 
ture in Fig.[5Ka). Two important features exist: (i) small 
particles reach the bottom of the flow faster compared to 
large particles reaching the top; (ii) large particles appear 
to rise predominantly together (indicated by the band of 
low (/>). The first feature is easily explained by asymme¬ 
try: small particles beginning the experiment near the 
interface between the two species quickly travel to the 
bottom through the large particle matrix, in accordance 
with P{Az s \(f>). The second feature is possibly linked to 
the large particles having a maximum segregation speed 
at an intermediate concentration. 

Theoretical comparison. —Current approaches to mod¬ 
eling size segregation use an advection-diffusion equation 
for (j) [icj ]: 

^ + div(0u ) -±( gm )=±( D fj, (1) 

where u is the bulk velocity field, q = q(u>^o,g) is the 
mean segregation speed, g is gravity, D is the diffusivity 
and F(<j>) is the flux function, which determines the de¬ 
pendence of the segregation flux on (j). The simplest flux 
function has a quadratic form 


Fit) = <K1 - 0)- (2) 

This is employed in a number of models 12 ■E3 and is 
symmetric about <fi = 0.5, dictating that small and large 
particles behave identically, but in opposite directions. 
Recently, asymmetric flux functions were introduced by 
[4l| with the simplest being a cubic form 


F((l>) = A K (j>(l - 0)(1 - K<f>), 


(3) 
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FIG. 5: (color online), (a) Temporal development of 4>{z,t) 
versus normalized flow height z/h for a <!> = 50% mixture with 
dmax = ±30°. (b)-(c) Theoretical predictions from Eq. (|4)>. 
(b) Prediction using the symmetric flux function ((2j, with 
S r = 0.016 and S r /D r = 20.9 [39]. (c) Prediction using the 
asymmetric flux function J3]), with S r = 0.030, S r /D r = 29.6 
and k = 0.89. 


where asymmetry parameter 0 ^ n < 1, and normaliza¬ 
tion constant A K gives the same amplitude as the sym¬ 
metric flux function. The applied shear gives a velocity 
profile u = ( u(z , t ), 0,0). In combination with the lateral 
spatial uniformity of this means that the transport 
term in Eq. m is zero. Equation m reduces to: 


d(f) 

~dt 


^(SrF^) 


— (D W) 

dz ^ r dz ): 


( 4 ) 


where z = z/A , and S r = qT/A , D, = DT/A 2 are 
non-dimensional segregation and diffusive-remixing co¬ 
efficients respectively. The symmetric and asymmetric 
models were least squares fitted to the data in Fig.[2][b) to 
obtain S r = 0.016 for the symmetric model and k = 0.89, 
S r = 0.030 for the asymmetric model. Integrating Eq. Q 
gives the cj) evolution in Figs.[5](b) and[5](c). Qualitatively, 
Fig.EIc) reproduces the experimental result on some crit¬ 
ical points: (i) the difference in time between the arrival 
of small particles at the bottom and large particles at the 
top of the flow; (ii) the collective rising of large particles; 
and (iii) a lower cf> in the bottom half of the flow near the 
end of the experiment, indicating that some large parti¬ 
cles are still inside the small particle matrix, segregating 
very slowly. These features are not found in the symmet¬ 
ric result in Fig.[5](b). The theoretical displacements per 
cycle are given by 


A zi = S r 


PM 

1 - 0 ’ 



( 5 ) 


and are shown alongside the experimental data in 
Figs. [4](c) and [U(d). The trend is clearly better pre¬ 
dicted by the asymmetric flux, which is able to reproduce 
both the peak in |(A^;)| around <f> = 0.5 and the nonlin¬ 
ear decrease of |(Az s }|. We attribute the discrepancy of 


|(Az s }| at low cf> to tracking errors, when small particles 
move more than their radius and their displacement is 
not recorded, thereby lowering the measured value. 

Discussion .—We analyze particle motion in a segre¬ 
gating bidisperse mixture under oscillatory shear and dis¬ 
cover an underlying asymmetry in the behavior of large 
and small particles. The small particle motion is step¬ 
like, falling down through the large particle matrix typ¬ 
ically one layer at a time. On the other hand, the large 
particle motion is smoother but slower, and linked to the 
scale of rearrangements of the surrounding small parti¬ 
cles. The asymmetric motion of the large and small par¬ 
ticles combine to give a characteristic dependence of the 
particle segregation speeds on the local volume fraction. 
Large particles segregate slower in the presence of many 
small particles, while small particles segregate faster in 
the presence of many large particles. We also observe 
that large particles move quickest when close to other 
large particles at intermediate concentrations, a process 
reminiscent of a collective motion (4^ |. The underlying 
asymmetry also manifests at rneso and bulk scales. In 
the development of </>(z,f), small particles reach the bot¬ 
tom of the flow faster than large particles reach the top. 
The segregation time increases linearly when a mixture 
contains a larger fraction of small particles. Although 
there is no direct evidence that the observed asymmetry 
persists for continuous shear, Staron and Phillips [ 20 J re¬ 
port that the segregation time under steady shear also 
increases linearly with the total concentration of small 
particles. These insights give a new understanding of 
segregation in sheared systems, with the dynamic behav¬ 
ior of two species being inherently different. 

Models for segregation have typically considered the 
motion of the large and small grains to be identical. 
However, an experimentally determined cubic flux 41] 
brings asymmetric behavior for the two species and gives 
good agreement on both particle and bulk scales. This 
draws parallels with the use of asymmetric flux func¬ 
tions to model asymmetry in sedimentation [43|, traffic 
flows 44], |45| and diffusion across membranes [46[ . For 
example, in the sedimentation of suspensions, particles 
settle faster when traveling together, but the settling ve¬ 
locity goes to zero more rapidly than a linear decrease at 
high concentrations [43]. Similarly, the velocity of cars 
in traffic approaches zero non-linearly at high car den¬ 
sities 0]. The commonality between these processes is 
their discrete nature, but interestingly, size segregation 
is the only process that consists of two discrete species. 

Further work is needed to analyze particle scale mo¬ 
tion for continuously sheared flows, e.g. down chutes 
and within rotating drums, to determine whether asym¬ 
metry persists and what form the segregation flux takes. 
The distinct segregation dynamics of the two species also 
leads to questions as to a possible relation with other dy¬ 
namic processes such as dynamic heterogeneity 48]. 
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